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We report the detection and observed characteristics of giant pulses from the 

pQ I Crab Nebula pulsar (B0531+21) in four frequency bands covering 20-84 MHz 

pn I using the recently-completed Long Wavelength Array Station 1 (LWAl) radio 

^ \ telescope. In 10 hours of observations distributed over a 72 day period in Fall 

of 2012, 33 giant pulses having peak flux densities between 400 Jy and 2000 Jy 
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^ I were detected. Twenty-two of these pulses were detected simultaneously in chan- 

ts ■ nels of 16 MHz bandwidth centered at 44 MHz, 60 MHz, and 76 MHz, includ- 

ed ' 

ing one pulse which was also detected in a channel centered at 28 MHz. We 

quantify statistics of pulse amplitude and pulse shape characteristics, including 
pulse broadening. Amplitude statistics are consistent with expectations based on 
extrapolations from previous work at higher and lower frequencies. Pulse broad- 

Q I ening is found to be relatively high, but not significantly greater than expected. 

^ I We present procedures that have been found to be effective for observing giant 

O ■ pulses in this frequency range. 
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Introduction 



The Crab Nebula pulsar (B0531+21) intermittently produces single p ulses having inten- 
sity o rders of magnitude greater than those of the normal periodic emission ( IStaelin &: Reifenstein 
19681 ). These giant pulses (GPs) consist of cornplex superpos i tions of nanosecond- and 
microsecond-scale impulses ( Hankins et al.l l2003l : ICrosslev et al. |2010[) and have been ob 



served over many octaves of the ra dio spectrum (ISallmen et al.l 1 19991 : ICordes et al.l |2004J : 



Popov et al.l I2OO6I : iBhat et al.l 120071 ). The details of the GP emission mechanism remain 
mysterious. Because the intrinsic duration of Crab GPs is much less than the pulse broad- 
ening due to the host nebula and the interstellar medium (ISM) along the line of sight, Crab 
GPs are also potentially of great value as a probe of the electron density of those regions. 

Studies of Crab GPs are often limited to frequencies above 100 MHz where the emis- 
sio n is relatively stron g and/ o r existing large rad io telescopes a r e ava il able; see for exara - 
plelMoffett &: HankinsI Jl996h: ICordes et al.rj2004h : IPodov et aP J2006I): Ishat et al.l ( 2007h: 
Maiid et al.lfcoilh i lKaruppusamv. Stappers fc Le3j2012[ ): and lMickahger et al.N2012l ). Stud- 
ies of Crab GPs below 100 MHz have historically been limited by the lack of suitable large 
telescopes and the difficulties imposed by strong and diverse radio frequency interference 
(R FI) at these frequenci es. Observations of Crab GPs have been reported at 4 and 60 MHz 
by iKuzmin et al.l (|2002[ ) and more recently at 23 MHz by iPopov et al.l (120061 ) . A search of 
14 hours of data collected in a 2.5 MHz bandwidth around 38 MHz by the 10-dipole Eight 
meter- wavelength Transient Array (ETA) experiment revealed 11 GP candidates, but these 
proved difficult to confirm due to the low sensitivity (~ 4cr for the strongest pulse candidate 
of ~ 1 kJy) and difficulties with RFI f lDeshpandell2009h . Detection of Crab GPs by LOFAR 
in the frequency r anges 32-80 MHz an d 139-187 MHz has been reported but further details 
are not provided (IStappers et al.l 1201 ll ). 



In this paper we report the detection and characterization of Crab GPs in the frequency 
range 20-84 MHz by the recently-completed Long Wavelength Array Station 1 (LWAl) radio 
telescope. LWAl consists of a pseudorandom array of 256 dual-polarized broadband dipoles 
within a 110 m x 100 m elliptical footprintjj and is capable of Nyquist-rate acquisition of 
beam outputs covering its entire tuning range of 10-88 MHz simultaneously. This facilitates 
high-signal-to-noise ratio detection of Crab GPs at a rate of several per hour; also, the large 
instantaneous bandwidth provides the ability to study variations in flux density and pulse 
broadening with frequency on a pulse-by-pulse basis. The large bandwidth also facilitates 
improved methods for the recognition and mitigation of the effects of RFI. 



^Not including a small number of "outrigger" dipoles which lie > 300 m from the center of the array and 
are used for array calibration and power spectrum measurement. 
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This paper is organized as follows: Section |2] describes LWAl and the methods used for 
data collection, data reduction, identification of Crab GPs, pulse characterization, and flux 
density calibration. Section [3] summarizes the findings of the study including statistics of 
pulse amplitude and pulse broadening. Section H] includes a brief description of ongoing and 
planned follow-up work. 



2. Instrumentation &; Data Analysis 

2.1. LWAl 

LWAl is collocated with the Very Large Array (VLA; 107.63° W, 34.07° N) in cen- 
tral New Mexico. The telescope consists of a primary array of 256 pairs of dipole-type 
antennas whose outputs are individually digitized and formed into beams. A detailed tech- 
nical description of th e instrument, including performance characterization, can be found in 



EUingson et al.l ( l2013l ). LWAl provides 4 independently-steerable dual-polarization beams. 
Each beam provides two separate passbands with independently-selectable center frequen- 
cies. Each passband consists of Nyquist-rate time-domain output at sample rates up 19.6 
million samples per second (MSPS) per passband with ~ 17 MHz usable bandwidth, in 
4-bits real ("I") plus 4-bits imaginary ("Q") format. The beam main lobe width varies with 
frequency and pointing direction, but ranges from about 23° full width at half-maximum 
(FWHM) at 20 MHz to about 3° FWHM at 84 MHz for observations of the Crab Pulsar 
within about 1 hour of upper culmination. The system equivalent flux density (SEED) varies 
as a function of frequency, pointing relative to zenith, and (because the system temperature 
is strongly dominated by Galactic noise) pointing in celestial coordinateso SEED estimates 
for the flux density calibration specific to the present work are developed in Section 12. 6j 
however SEED is typically found to be in the range 6-17 kJy for high-elevation pointings 
over most of the frequency rangejj This corresponds to a 5a sensitivity < 100 Jy assuming 
both polarizations, 17 MHz bandwidth, and 25 ms integration. 

A descr iption of LWAl ear ly science results (including first pulsar observations) can 



be found in iTaylor et al.l ( 120121 ). To date LWAl has detected 13 known pulsars and has 



several on-going programs dedicated to survey, monitoring, and characterization of pulsars 



^SEFD is defined here as the flux density of an unresolved source in the main lobe which results in a 
doubling of the measured beam output power relative to the value in the absence of the source. 

■^The counterintuitive fact that SEFD is approximately constant over a ~ 4 : 1 frequency range can be 
understood by noting that the collecting area of a dipole decreases as the square of frequency, whereas the 
Galactic noise-dominated antenna temperature decreases with frequency at about the same rate. 
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(jStovall et alJl2013l ) 



2.2. Data Collection 

Data for the present work were collected in 10 1-hr observations in Fall 2012, each 
centered on the transit of the Crab Pulsar, as detailed in Table [H In each observation, two 
beams were used together to track the pulsar (i.e., both beams pointing in precisely the 
same direction). The center frequencies of the passbands of the first beam were 76 MHz and 
60 MHz, whereas the center frequencies for the second beam were 44 MHz and 28 MHz. Since 
the usable bandwidth is > 16 MHz per passband, this scheme yields contiguous frequency 
coverage from 20 MHz to 84 MHz. The total volume of data collected in this fashion 
was ~ 5 TB; relatively large since the data are recorded as raw voltage streams without 
channelization or other "on the fly" processing. This format provides the maximum flexibility 
for both initial processing and subsequent reprocessing of the data at later time using other 
search parameters or techniques. 

Although the instrument is capable of tracking sources to very low elevation angles, 
we have restricted the observations to 1 hour each day around upper culmination in or- 
der to (1) limit the extent to which change in the beam characteristics over the track can 
manifest as variations in estimated flux density, (2) distribute the observations more evenly 
over the expected range of source conditions, and (3) reduce the tempo of data acquisi- 
tion to a manageable pace. Concerning (2), a parameter of interest in the present work 
that is known to vary s ignificantly over timescales of weeks to months is pulse broaden- 



i ng (IKuzmin et al.l 120081); in fact, such variat ions can occur even on timescales of minutes 



(jKaruppusamy. Stappers fc van StratenI 120101 ). 



Because this study involved the use of a new instrument and new data analysis pro- 
cedures, we performed simultaneous observations of the bright pulsar B0329+54 using the 
remaining two beams, using the same synchronized-tracking four-passband scheme employed 
for Crab GP observing. This allowed an independent check of the instrument and data reduc- 
tion procedures using a well-known object, as well as providing an "off pointing" (separated 
by ~ 2 h in right ascension) for confirmation of Crab GP detections as well as for identifica- 
tion of RFI through anticoincidence techniques. 
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2.3. Data Reduction 

The data were analyzed using custom incolierent dedispersion software developed by one 
of us (S.E.), designed specifically for rapid analysis of periodic and single-pulse emission of 
dispersed astrophysical signals in raw LWAl beam data. The steps in data reduction for one 
passband of one beam were as follows: First, any samples having a//^ + Q^ > 7 (that is, the 
magnitude of the complex- valued sample equal to or greater than maximum encodable real 
or imaginary magnitude) were assumed to be due to RFI and set to zero; typically 0.1% to 
1% of samples were zeroed in this process. Next, data were partitioned into blocks of length 
L (center frequency- dependent) and divided into channels using the fast Fourier transform 
(FFT). Values of L and the associated time and frequency resolutions are shown in Table [2J 
These values were chosen to make the residual dispersion negligible relative to the width of 
the scatter-broadened pulse after incoherent dedispersion. The raw linear polarizations were 
incoherently combined at this point. Time-frequency pixels associated with the outer 10% of 
frequency channels were discarded and the inner 90% were "flattened"; that is, variations in 
frequency and time (attributable to the instrument or to ionospheric effects) within the time 
span of the spectrogram (typically 10-50 s) were removed. Finally, all pixels in frequency 
bands containing persistent RFI, or RFI which was intermittent but frequently observed, 
were discarded; additional details about this procedure are given in Section 12.41 

In addition to the amplitude variation, the ionosphere imparts both refraction and 
dispersion to received signals. For the purposes of the present work, these additional effects 
are negligible and may be safely ignored: T he worst-case ionospheric refraction is on the order 



of arcminutes (see e.g. iKassim et al.l (120071 )) and thus is a tiny fraction of the narrowest beam 



width in the present work. The worst-case ionospheric contribution to the total dispersive 



delay over 20-84 MHz is ^ 1 ms (JRoussel-Dupre. Jacobson fc Triplettll200ll ). corresponding 



to an apparent error in dispersion measure (DM) < 10~^ pc cm~^. This is at least two orders 
of magnitude less than the error in assumed DM used for dedispersion, which is described 
below and shown to be insignificant. 

Following R FI mitigation, the resu lt was incoherently dedispersed assuming DM = 
56.791 pc cm~^ (JManchester et al.l 120051 ). Although coherent dedispersion would be prefer- 



able from a sensitivity perspective, the residual dispersion after application of incoherent 
dedispersion at the time-frequency resolutions identified above is less than than the ob- 
served width of the scatter-broadened pulses by more than two orders of magnitude; thus 
the penalty in the present work is negligible. In addition to having a large advantage in 
terms of computation burden, the frequency- domain interference mitigation procedure de- 
scribed below (Section 12.41) is easily implemented in the incoherent dedispersion framework, 
but much less straightforward to implement as part of a coherent dedispersion scheme. The 
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dedispersed time series was averaged to obtain an integration time and effective time reso- 
lution of approximately 25 ms for the 76 MHz and 60 MHz passbands, and either 33 ms or 
50 ms for the 44 MHz and 28 MHz passbands. 

Note that there is no subsequent adjustment to the assumed DM of 56.791 pc cm~^ 
in our processing. In particular, we did not perform a search over DM. As explained in 
Section 13.21 refinement of the DM used for dedispersion proves to be unnecessary because 
the resulting bias in estimated pulse parameters turns out to be negligible. 

The products saved from this procedure were (1) the dedispersed and averaged time 
series, and (2) the associated spectrogram, which was further averaged to a time-frequency 
resolution of approximately 5 kHz x 100 ms. A representative spectrogram is shown in 
Figure [U The spectrograms play an important role in data interpretation, as explained 
below. 

In observing sessions where we simultaneously observed B0329+54, we repeated the 
above procedure in all details, except for the DM (26.833 pc cm~^) and we t ook the additional 



step of folding the resulting time series at the known period (~ 0.714 s) (JManchester et al. 



20051 ). In each case we detected the pulsar in all four passbands with high signal-to- noise 
ratio, giving confidence in the new instrumentation and data analysis procedures. Further- 
more we confirmed that GPs detected in the beams pointing to B0531-I-21 were absent from 
the beams pointing to B0329+54, as expected. 



2.4. Identification of Crab GPs 

Traditional methods for identification of GPs involve matched filtering of the time series, 
typically with iteration over DM in order to confirm that detected pulses are in fact Crab 
GPs. A straightforward application of this procedure is not robust for the present work, for 
several reasons. First, ionospheric scintillation results in variations that can vary by as much 
as 10% of the system temperature with characteristic time scales on the order of minuteso 
This frequently results in large variations in the noise baseline of the time series, and occurs 
on time scales roughly equal to or less than the time it takes for the "chirp" waveform of 
the dispersed pulse to traverse the passband. Furthermore, GPs in this frequency range 
remain scatter-broadened to durations on scales of 100-1000 ms after dedispersion. This 
tends to make matched filter detection procedures that assume canonical pulse shapes with 
time-invariant noise baselines unreliable and frustrating to use, especially in the presence of 



*For examples, see Figures 10-11 of lEllingson et al.l (|2013l ). 
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impulsive RFI. 

Furthermore, impulsive RFI is ubiquitous in this frequency range. In particular, we have 
found that signals from terrestrial television transmitters in the frequency range 54-87 MHz 
are a serious problem. Transmissions in both the modern digital "ATSC" format currently 
in widespread use in the U.S., as well as transmissions in the legacy analog "NTSC" format, 
are routinely found in LWAl data. ATSC signals are recognizable as emissions having 6 MHz 
instantaneous bandwidth, whereas NTSC signals are recognizable as relatively narrowband 
video and audio carriers separated by 4.5 MHz. With the exception of a persistent NTSC 
video carrier seen at 55.25 MHz (associated with a station in Mexico), these emissions 
are usually not seen as continuous transmissions, but rather in bursts with durations on 
the order of milliseconds to seconds. These bursts are believed to be due to reflections 
from ionized meteor trails and low Earth orbiting satellites, as opposed to reflection from 
aircraft or the ionosphere|f| These emissions have two particularly irritating features from 
the perspective of searching for single dispersed pulses: (1) They frequently exhibit the same 
"reverse ramp" waveform expected of scatter-broadened pulses, and with similar durations, 
and (2) Because they are narrowband relative to the overall bandwidth, they are often not 
significantly suppressed by the process of dedispersion - they are simply translated in time 
and (in the case of ATSC flares) also "smeared" in time. 

Excluding bandwidths of 1.0 and 0.5 MHz around frequencies of NTSC video and audio 
carriers (respectively) dramatically improves - but does not eliminate - the NTSC problem. 
Specifically, strong NTSC flares occasionally "overrun" the excluded frequency ranges. ATSC 
flares cannot be effectively mitigated by frequency blanking without removing the entire 
6 MHz bandwidth of the emission. Reducing the effect of NTSC and ATSC flares to a 
negligible level by frequency-domain blanking requires eliminating an unacceptable amount 
of otherwise-usable bandwidth. Because the flares occur with relatively small overall duty 
cycle, an automated joint time-frequency blanking scheme probably would be acceptable; 
however, this requires significant software development and has not yet been implemented. 

At present, we find that the most practical and robust method for detection of Crab GPs 
in our frequency range is as follows: First, we manually search the spectrograms for each pass- 
band for the characteristic "chirp" of a Crab GP which is visible over most of the passband. 
Figure [1] shows an example. We then confirm the presence of the distinctive broadened pulse 
waveform at the appropriate time in the dedispersed time series. Typically - but not always 
- pulses are first detected in the 76 MHz or 60 MHz passbands; this is because the pulse 



'^This is based on a preliminary and ongoing analysis of LWAl data by J. Helmboldt of the U.S. Naval 
Research Laboratory. 



broadening is less severe, and becau se pulse flux density in this range of frequencies generall y 



increases with increasing frequency (jPopov et al.ll2006l : iKaruppusamy. Stappers fc Leell2012l ) . 



The dedispersed time series for the other passbands is then also checked, accounting for the 
expected dispersion delay which is tens to hundreds of seconds due to the large fractional 
bandwidth between passbands. The first 10 or so GPs that we detected using this approach 
were subsequently independently confirmed by repeating the procedure while sweeping in 
DM. We soon realized that that this computationally-expensive procedure was not neces- 
sary for the present work, since simpler procedure of detecting pulses through spectrograms 
and/or appearance in multiple passbands with the expected dispersion delay yielded a false 
alarm rate of zero, whereas the results of DM sweeps were frequently difficult to interpret 
due to the effects of the ionosphere and impulsive RFI, as discussed above. 

A possible disadvantage of this detection procedure is somewhat reduced sensitivity. 
Manual inspection of dedispersed time series reveals many pulse candidates which are quite 
likely to be Crab GPs, but which are not visible in spectrograms and/or ambiguous or non- 
detections in the dedispersed time series for other passbands. We suspect refinements to 
the above procedures can increase the rate at which GPs are detected by a factor of 2-3 at 
least. For this initial work, however, we prefer a simpler and more conservative approach. 
Furthermore (as discussed in Section [23]) , we shall apply an even more stringent criterion 
of S/N > 4.4cr for the peak flux density. Several candidate GPs detected using the above 
procedure failed this criterion, indicating that the sensitivity penalty is actually not unusually 
severe in comp arison to se a rch p roced ures (including fully-automated techn iques) employed 



by others; e.g. iBhat et al.l (120071 ) and IKaruppusamy. Stappers fc Led (120 12[ ). 



Finally, we note that frequencies below 54 MHz down to the lowest frequency consid- 
ered in this study (~ 20 MHz) exhibit relatively sparse and low levels of anthropogenic 
RFI during local hours of darkness, when all of our observations were conducted. We found 
that no frequency blanking was required for the 28 MHz passband, and only 0.2 MHz was 
blanked in the 44 MHz passband (to remove a persistent tone - possibly locally generated 
- at 48.0 MHz). Ironically, ionospheric scintillation is significantly worse at the lower fre- 
quencies; in some cases the resulting variations in the noise baseline were large enough to 
make characterization of the detected pulses intractable, and so these pulses were excluded. 
Nevertheless it is believed that it is the low fiux density of pulses relative to the system sensi- 
tivity, as opposed to the disposal of pulses distorted by baseline variations, that is primarily 
responsible for the relatively low rate of detections in the 28 MHz and 44 MHz passbands. 
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2.5. Pulse Characterization 



Below 84 MHz, the duration of Crab GPs is broadened by at least three orders of 
magnitude by scattering within the host nebula and along the line-of-sight through the ISM. 
Thus it is an excellent approximation to assume that the shape of the received pulse is simply 
a s caled replica of th e scattering impulse response g{t). Appropriate pulse shapes are given 
by IWilliamsonI ( 1l972l ) as follows: 



git) 



oc 



'l/Td)exp{-t/Td)u(t) 



corresponding to scattering by a thin slab of material of infinite extent, and 

git) oc (7rV]/4t^) exp (-7rVrf/4t) M(t) 



1^ 



(2) 



corresponding to scattering by material which is uniformly-distributed along the line of sight. 
In both forms u{t) is the unit step function {u(t) = for t < 0, u(t) = 1 for t > 0; t = 
taken to be the pulse start time) and t^ is referred to as the characteristic broadening time. 
The dependence of Td on frequency in the first form is determined by the spatial distribution 
of electron density. When this distribution follows the Kolomogorov s pectrum, r^ oc t^~^'^, 
and e arly measuren i ents seemed to indicate that this was the case (jlsaacman &: Rankin 



19771 : ISallmen et al.l Il999l ). Subsequent observations indic ate a considerably more shal- 



low dependence, with exponents in the ra n ge —3.5 to —3.2 (iPopov et al.l 120061: iBhat et al. 



20071 : iKaruppusamy. Stappers &: Leell2012l ) . iKaruppusamy. Stappers fc Led ( 12012| ) note that 



pulses observed in the 116-174 MHz band are not well described by Equation [H but obtain 
good fits to the modified form 



g{t) (X t^exp{-t/Td)u{t) 



(3) 



which includes a rounded leading edge with rise time determined by the parameter /3. We 
also find in the present work that the above form provides the best overall fit to the observed 
pulses, whereas the other forni s gene rally do not; thus we have adopted this model. As in 
Karuppusamy. Stappers fc Led (|2012| ). we constrain /3 to be in the range 0-1 and find that 
this range seems to accommodate all pulses encountered. 

Thus, each pulse in the present work is characterized by r^, /3, and amplitude. These 
are determined as follows: First, we fit a second-order polynomial to the time preceding 
the start of the pulse (i.e., t < 0) in the dedispersed time series. The baseline variation is 
removed, and then an automated brute-force fit in three dimensions {ra, /3, and amplitude) 
is performed to find the least-squares fit to the data. Figure [2] shows examples. 

For verification, the model pulse is then subtracted from the data to confirm that the 
resulting time series is fiat and free from artifacts. This turns out to be particularly important 
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in our frequency range, since detectable pulses sometimes overlap each other and thus bias the 
estimation of parameters. In the present work, two otherwise-qualified pulses were excluded 
from consideration (and thus do not appear in Table [T]) for this reason. 

The amplitude of the pulse is characterized in two ways. First, the pulse peak S/N is 
calculated as the ratio of the peak value of the model pulse to the pre-pulse noise power 
A^. This is the value used in the 4.4o" detection criterion, described above. Because the 
integration time At (25-50 ms) is at least an order of magnitude less than the apparent 
width of the pulse, the effective S/N is much greater, and is given by 



[{EJW) /N] ^M[W,. 



int 



(4) 



where Eg is the pulse energy, W is t he "effective widt h" of the pulse, and Wmt is the intrinsic 
pulse width, assumed to be 300 /xs (JBhat et al.ll2007l ). The effective width W is the width of 
a rectangular pulse having the same peak amplitude and energy as the model pulse. For the 
form given in Equation [3l W ^ Td exp {\fP). Note that the effective S/N can only be realized 
for detection by eithe r optimal pulse-matched filtering or by deconvolution of g{t) from the 
intrinsic pulse shape (jBhat. Cordes &: Chatterjedl2003l ). Thus the estimated value obtained 
using the above procedure is used only for flux density calibration (and not detection), as 
described in the next section. 



2.6. Flux Density Calibration 

Flux density calibration is challenging for array radio telescopes operating at low fre- 
quencies. The primary difficulty is that beam characteristics vary as a function of frequency 
and pointing relative to zenith, pointing relative to the celestial sphere (because the system 
temperature may be dominated by sky noise), an d source magnitudes ca n be modulated by 
as much as 10% due ionospheric scintillation. See lEUingson et al.l ( 120131 ) for a discussion of 
these factors as they apply to LWAl. 

For the present work we used the tr a nsit d rift scan SEFD measurement procedure 
described in Appendix B of lEUingson et al.l (120131 ). In this procedure, the data are collected 
in exactly the same manner described in Section 12.21 with the exception that the beams are 
fixed, pointing at the point of upper culmination of the pulsar, and the total duration of the 
observation is 2 hours. The data are reduced to full-bandwidth total power (single channel) 
time series. The transit of Tau A (the unresolved radio source which includes the Crab 
Pulsar and the surrounding nebula) is clearly visible as a broad peak in the time series as 
the source moves through the main lobe. When the source is outside the main lobe, the time 
series is thus a representation of the system temperature of the instrument in the vicinity 
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of, but excluding Tau A. The peak of the time series can thus be interpreted as the sum of 
the SEFD and the flux density of Tau A. Assuming the flux density of Tau A is known, the 
SEFD can in principle be calculated from the ratio of the peak of the time series to the level 
of the flat, off-peak portion of the time series. In practice, we estimate and subtract a model 
of the beamshape, and use the level of the resulting noise baseline at the time of transit as 
the "source turned off" value; this reduces the error associated with the gradual change in 
SEFD over time. 

In the present work we ignore the variation in beam characteristics with source position 
(in particular, source elevation). Since each observation is limited to 1 hour around upper 
culmination, this is variation is believed to be negligible. For example: Assuming a cos^ Z 
dependence of beam directivity with zenith angle Z, the associated change of directivity over 
a 1 hour observation is less than 2%. 

The flux calibration procedure requires an accurate value for the flux density of Tau A 
over the range of frequencies considered. Our model for the flux density of Tau A is 



(1944 Jy) (Z//76 MHz) 



-0.27 



(5) 



The spectral index in this model is the one derived in lBaars et al.l (119771 ). The reference flux 
density of 1944 Jy at 76 M Hz is the average of values ext rapo lated from m easurements at 
22.25 MHz and 81.5 MHz bv lRoger. Costain fc LacevI Jl969h and lParkeil Jl968h . respectively0 
using this spectral index. This model gives values which fall within the error bars of the 
22.25 MHz and 81.5 MHz measurements (±6% and ±4%, respectively) and turns out to 
be essential l y the s ame as model used in recent wo r k on C rab GPs at low frequencies (e.g., 



Bhat et al.l (120071 ): iKaruppusamy. Stappers fc Led (120 12l )). attributed to iBietenholz et al. 



(119971 ). which uses the same spectral index but gives a 76 MHz flux density of 1915 Jy. 
Uncertainty in the reference flux density (±5%) and the spectral index (±0.04), together 
result in an uncertainty of about ±10% at 20 MHz relative to the 84 MHz value. In practice 
the model may also be biased at the lowest frequencies by factors not taken into account; 
e.g., turnover in the spectrum of the nebula. 

Using the above approach we obtain SEFD of 13.7 kJy, 10.6 kJy, 10.3 kJy, and 11.5 kJy 
for the 76, 60, 44, and 28 MHz pass bands, respectively. T he 76 MHz value is approximately 
twice the 74 MHz value reported in lEUingson et al.l (120131 ) . which was obtained using essen- 
tially the same technique. The difference is due to a recalibration of the array beamforming 
delays which occurred before the observations reported here were started, and appears to 



^These measurements are also reported in lBaars et al.l (|l977l ) 
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have degraded the beamforming performancei] No earlier measurements of the same point- 
ing are available for the other passbands; however the values reported above fall within the 
range 6-17 kJy that we normally expect for high-elevation sources. The clear presence of a 
peak at the expected time in the flux density calibration drift scans confirms pointing was 
working properly for all passbands during the observations. 

From the above SEFD estimates we find the 4.4cr sensitivities (in terms of the effective 
S/N as described in the previous section, and not in terms of the peak waveform S/N) are 656, 
515, 444, and 493 Jy for the 76, 60, 44, and 28 MHz passbands, respectively. These figures 
take into account frequency- domain blanking for RFI purposes, resulting in bandwidths of 
14.07, 13.64, 17.44, and 17.64 MHz for the 76, 60, 44, and 28 MHz passbands, respectively. 



3. Results 

The observing sessions and the resulting detections are summarized in Table [T] A total 
of 33 Crab GPs were identified in 10 hours of observation, distributed over 72 days during 
Fall 2012. One pulse was detected simultaneously in all 4 passbands, 21 pulses were detected 
simultaneously in 3 passbands, and an additional 11 pulses were detected simultaneously in 
2 passbands. All but a handful of pulses detected passed the 4.4cr peak S/N criterion. 



3.1. Pulse Amplitudes 

Traditionally, statistical distributions of pulse amplitudes are expressed as power law 
fits to the tail distributions of th e corresponding p r obabi li ty distribution fuii ctions (PDFs, 



or "histograms" ) ; see for example lArgvle fc Gowerl ( 119721 ) ; iMajid et al.l (l201l[ ) ; and 



Karuppusamy. Stappers fc Led ( 120121 ) . This form of analysis is not meaningful for the present 
work due to the relatively small number of pulses detected. Instead, Figure El summarizes 
detections as amplitude cumulative distribution functions (CDFs); i.e., expected rates of de- 
tections for pulses > a specified flux density. The CDF is the integral of the PDF; however 
the tail distributions of the two functions are not related in a straightforward way. Never- 
theless for the purposes of pulse rate prediction in this and future work, a power law fit to 
CDFs shown in Figure [3] was done, with the following results: 

N,, = (1.08 h-i) {SMjY^-'' (6) 



^A recalibration of the array in March 2013 has restored the original sensitivity. 
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iVeo = (0.47 h-^) (5/kJy)-^-^^ (7) 

iV44 = (0.25 h"^) {S/kJjy^-^^ (8) 

Where N^, is the number of pulses expected with flux density > S a.t u MHz. (No such 
characterization is possible for the 28 MHz passband, since one pulse was detected.) Note 
that the distributions exhibit very similar dependence on frequency, and may be steepening 
with increasing frequency. Although these values cannot be directly compared to analogous 
characterizations of the PDFs, we note that the same trend is seen in past work using the 
PDFs. 

Figure H] shows the brightest pulses detected in each passband, the system sensitivity to 
Crab GPs (worked out at the end of Section [2^ . and predicted flux densities for detection 
rates of 1 pulse/hour using Equations [UHHl Because the brightest pulses represent a under- 
sampled ta i l distr ibution, it is not reasonable to infer a spectral index from our data alone. 
Bhat et al.l (120071 ) suggest a power law of ~ z/+^-^ for flux de nsity from a two-point fit using 



the 23 and 200 MHz data. Although the 174 MHz value from lKaruppusamy. Stappers fc Lee 



( I2OI2I ) is for only 1 h, a model for the tail distribution of the PDF is also available; this 
model predicts that the brightest pulse in 10 h (for comparison with the present work) is 
10.5 kJy. This extrapolation combined with our data suggest that the dependence of flux 
density on frequency may be significantly less than z/^^'^. 



3.2. Pulse Broadening 

Table |3]shows statistics of r^ and the pulse rise time parameter /3 for all pulses considered 
in the study. Figure \5\ shows the statistics of r^ in context with the analogous values for 
previous studies at other frequencies, which are also listed in Table El As explained in 
Section 12.51 the exponent a of the power law dependence of r^ with frequency is of particular 
interest as it is related to the distribution of the scattering material. Table H] shows a 
calculated using various combinations of the data represented in Figure [5] and Table |3l Data 
from our study seem to indicate that (1) t^i is significantly higher than might be expected 
by extrapolation from other studies, and (2) the dependence on frequency is significantly 
weaker than might be expected by extrapolation from other studies. 

With respect to the first finding, we have ruled out errors in assumed DM and residual 
dispersion (due to the use of incoherent dedispersion) ; we find that no combination of these 
mechanisms can account for a bias greater than about 35 ms in r^, which is less than 8% of 
the observed mean value for the highest-frequency passband. The largest contribution to this 
bias is 28 ms due to error in the assumed DM. This error is quantified as follows: The Jodrell 
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Bank Crab Pulsar ephemeriqfl ( Lyne. Pritchard fc Grahani-Sniithlll993l ) indicates that the 
actual DM for the months of October 2012 through January 2013 was consistently higher 
than the assumed DM, with the worst case error being 0.0644 pc cm~'^ larger than the 
assumed DM. This corresponds to an error in dispersive delay of about 636 ms over the 
entire frequency range 20-82 MHz. However we process this frequency span in separate 16 
MHz passbands, and the worst case delay error across a passband is only about 109 ms (in 
the 44 MHz passband; it is worse in the 28 MHz passband but no pulse width measurements 
were obtained in that passband). The bias subsequently induced in the estimate of r^ is only 
about one-fourth of 109 ms (i.e., 28 ms), which was determined as follows: We performed a 
simulation in which a pulse following the model of Equation |3] with /3 = 0.71 and Tfi=978 ms 
(mean values for pulses detected in the 44 MHz passband) was "smeared" by an error of 
109 ms in dispersive delay, and then we used our model-fitting algorithm to determine the 
apparent values of (3 and r^. We find that the apparent value of (3 (as well as pulse amplitude) 
is not significantly affected and that the apparent value of r^ is high by about 28 ms. The 
excess delay and resulting bias in all passbands is shown in Table |2l 

It is more likely that our data reflect actual changes in the distribution of ionized 
material within the nebula: The Crab Pulsar is known to exhibit increases in r^ by factors of 
2-10 for intervals on the order of days to mon ths, attributable to ionized cl o uds or filaments 
within the nebula crossing the line of sight ( iRankin fc CounselmanI Il973l : iLyne fc Thorne 



1975; Sallmen et al. 1999; Kuzmin et al. 2008 



Stappers et al.ll201ll ). This is consistent with 



the relatively high DM (noted above) in effect over the period of the observations. 

Power laws in the form r^ oc z/" implied by our data are included in Table HI Recent 
studies (summarized in this table and in Table |3]) have suggested that a lies in the range 
—3.5 to —3.2. This is not necessarily inconsistent with our data if we consider only the 76 
and 60 MHz data. However our 44 MHz data imply a flattening in the dependence with 
decreasing frequency which does not appear to be consistent with a < —3.2. The weak 
dependence on r^ on frequency implied by our data, like the large values of r^ implied by 
our data, might be interpreted as a perturbation du e to in homogeneity in the distribution of 
ionized material within the nebula. ICordes &: Laziol (l200l[ ) have shown that inhomogeneities 
in the form of disks and filaments may lead to large variations in the dependence of t^ on 
frequency, and that these effects may be particularly strong at low radio frequencies. 

Given that the studies identified in Table |3] and Figure |5] represent widely-separated 
epochs of scattering behavior within the nebula, it may be too much to expect a simple 
dependence of r^ on frequency to emerge from these data. Thus a broader interpretation 



*http://www.jb. man. ac.uk/^piilsar/crab. html 
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of our findings may be simply that the frequency dependence is complex: neither static 
nor well-described by a single power law extending to low radio frequencies. In any event, 
our results, like those of other recent studies, strongly suggest that the electron density 
distribution is not in the form of a Kolomogorov spectrum (a = —4.4), and in fact a (to the 
extent that it is reasonable to characterize the frequency dependence of pulse broadening in 
this manner) is probably always significantly greater than —4. 



4. Concluding Remarks 

We have demonstrated the ability of LWAl to observe single Crab GPs over a frequency 
range of ~4:1 at a rate of ~ 3 h~^, yielding information about pulse amplitude and pulse 
broadening that may be useful in gaining insight into the nature of the Crab GP emission 
mechanism as well as providing additional information about the electron density distribution 
within the nebula and along the line of sight. The key to continued progress is observing 
a far greater number of pulses distributed over a timeframe significantly greater than the 
weeks-to-months-long scattering epochs. To this end, we are continuing our program of 
Crab GP observations with the goal of increasing the number of detections by an order 
of magnitude and then repeating the analysis presented above. Also, we are preparing for 
simultaneous multiband observations of Crab GPs using LWAl with instruments operating 
at higher frequencies, which may provide additional information on Crab GP behavior by 
allowing us to simultaneously observe emission on both sides of the turnover in the GP 
spectrum. 

The authors acknowledge helpful discussions with W.A. Coles, T.H. Hankins, J.F. Helm- 
boldt, N.E. Kassim, W.A. Majid, and B.J. Rickett. Some of the data presented in this paper 
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University of New Mexico, Virginia Tech, and the Jet Propulsion Laboratory and in which 
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56291 
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56292 





2 


4 


3 


1 


3 


TOTAL 


1 


28 


33 


31 


22 


11 


> 4.4a 


1 


21 


32 


29 







Table 1: Summary of daily observing sessions. "MJD" is Modified Julian Date. Each 
session was 1 hour in length. Gaps of greater than 1 week are indicated by horizontal lines. 
The center four columns are the number of pulses detected in each passband. The rightmost 
two columns indicate the number of pulses appearing simultaneously in multiple passbands, 
respectively. The bottom row indicates the totals after excluding pulses having apparent 
peak S/N < 4.4a (25 ms, ^ 17 MHz). 



Passband 


Span 


L 


Az/ 


Ar 


Excess Delay 


Bias in r^ 


76 MHz 


68-84 MHz 


16384 


1196 Hz 


836 /is 


+20 ms 


+5 ms 


60 MHz 


52-68 MHz 


16384 


1196 Hz 


836 /is 


+41 ms 


+10 ms 


44 MHz 


36-52 MHz 


65536 


299 Hz 


3344 /is 


+109 ms 


+28 ms 


28 MHz 


20-36 MHz 


131072 


150 Hz 


6687 /is 


+466 ms 


+ 117 ms 



Table 2: Time-frequency resolutions used for incoherent dedispersion. Also shown is the 
additional dispersion delay across each passband due to worst-case error in the DM assumed 
for dedispersion (see Section 12. 3p , and the associated bias in r^ (see Section 13. 2p . 
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Fig. 1. — A pair of Crab GPs displayed in joint time-frequency ("spectrogram") represen- 
tation (top) and total power (bottom). The vertical axis of the spectrogram is a 19.6 MHz 
span centered at 76 MHz, with frequency increasing toward the top. Pulses begin at 2029 s 
and 2057 s at the top of the spectrogram, and "chirp" toward the bottom right. Horizontal 
white bands are excised frequency channels, including the outer 10% of the Nyquist band- 
width and (from bottom to top) NTSC channel 4 audio, NTSC channel 5 video, and NTSC 
channel 6 video. The vertical streaks at 2037 s and 2057 s are "flares" from ATSC channels 
6 and 4, respectively. The resolution of the spectrogram is approximately 100 kHz x 5 ms. 
The resolution of the time series is 836 /xs. 
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Fig. 2. — A pulse observed simultaneously in the 76, 60, 44, and 28 MHz passbands; top to 
bottom, respectively. Dispersion and dispersive delay between passbands has been removed, 
leaving only pulse broadening. The best-fit model pulses are also shown (note that most of 
the rising edge of the 28 MHz pulse is not visible due to inadequate sensitivity). The vertical 
scales are unmodified from the raw data, and thus include instrumental variations in gain 
and sensitivity between passbands. Date of observation: MJD 56291. 
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Fig. 3. — Detection rates expressed as number of pulses detected with flux density > the 
abscissa. Rates are expressed as the total number of detections divided by the number 
of hours of observation. The curves correspond to the (left to right) 44, 60, and 76 MHz 
passbands. (No data are shown for the 28 MHz passband since there was only one detection.) 
Also shown is a power law fit to the 18 largest values in each passband. 
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Fig. 4. — Red markers with error bars indicate the brightest pulses detected in the 28, 44, 60, 
and 76 MHz passbands in the present work. Also shown are the effective pulse sensitivities 
(4.4o") worked out in Section 12.61 and (for 44, 60, and 76 MHz only) flux densities at which 
a rate > 1 pulse h~^ is predicted, using Equations [6H8l Included f or compariso n are p oints 
representing the brightest pulse in 12 h at 23 MHz reported by Popov et al.l ( l2006l ). the 



at 200 MHz by 



Bhatetal 



fl2007f ) 



brightest pulse in 14 h at 38 MHz reported bv iDeshpandd (120091 ) . the brightest pulse in 
0.25 h at 1 11 MHz reported bv Ipodov et al.| J2006h . the brisrhtest pulse in 1 h at 174 MHz 



reported bv iKaruppusamv. Stappers fc Led (120121 ) . and the brightest pulse reported in 3.5 h 
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Fig. 5. — Characteristic pulse broadening time r^ for GPs detected in the 44, 60, and 76 MHz 
passbands in the present work is shown in red. The error bars indicate the mean and ±lcr 
values. Also shown are measurements from previous studies listed in Table [31 The lines are 
best-fit power laws using the mean values from all studies shown (upper line, a = —3.7) and 
from all studies excluding the present work (lower line, a = —3.5); see Table IH 
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V [MHz] 


Td [ms 


Reference 


Comment 


23 


3000 ± 1000 


Popov et al. (2006) 




40 


600 ± 250 


Kuzmin et al. ('2002) 




44 


978 ± 287 


This work 


/3 = 0.71 ±0.20 


60 


140 ± 70 


Kuzmin et al. f2002) 




60 


768 ± 273 


This work 


/3 = 0.56 ±0.16 


76 


439 ± 122 


This work 


/3 = 0.42 ± 0.24 


111 


15 ±3 


PoDov et al. f2006) 




115 


13 ±5 


Staelin & Sutton fl970) 




157 


3.8 ±1.3 


Staelin & Sutton (1970) 




174 


1.5 ±0.4 


Karuppusamy, Stappers & Lee (2012) 


See caption 


200 


0.670 ±0.100 


Bhat et al. f2007) 




300 


1.3 ±0.2 


Sallmen et al. (1999) 


See caption 


300 


0.28 


Sallmen et al. (1999) 


See caption 


406 


0.18 ±0.08 


Kuzmin et al. ('2002) 




594 


0.018 ±0.008 


Kuzmin et al. ('2002) 




600 


0.095 ±0.005 


Sallmen et al. (1999) 




600 


0.043 ± 0.005 


Popov et al. (2006) 





Table 3: Measurements of c harac teristic pulse broadening times shown in Figure |5l 
Karuppusamy. Stappers &: Led (l2012l ) reports values for 7 frequencies between 116.75 and 
173.85 exhibiting ^ i/~^-^ dependence. ISallmen et al.l ( 119991 ) reports 300 MHz values for two 
widely-separated scattering epochs. 



Data 



a 



All data -3.7 

All data except 76, 60, and 44 MHz values from this work —3.5 

KaruDDUsamv. Stappers fc Lee f2012) (117-174 MHz) -3.2 ± 0.1 

76 and 60 MHz values from this work (only) —2.4 ± 2.9 

76, 60, and 44 MHz values from this work (only) — 1.4 ± 1.1 



Table 4: Power laws for the dependence of r^ on frequency inferred from various combinations 
of the data in Table [31 Shown is the exponent a for r^ oc z/" as determined by log-linear 
least squares fit to the indicated data. 



